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INTRODUCTION 


Transonic cryogenic wind tunnels, such as the Langley 0.3-Meter 
Transonic Cryogenic Tunnel (TCT) and the National Transonic Facility 



to either heterogeneous (pre-existing seed particles) or homogeneous 
(the gas forms its own seed particles) nucleation limits the minimum 
operating temperature and, consequently, the maximum Reynol ds number 
capability. This minimum operating temperature can be conservatively 
limited to temperatures corresponding to saturati on at the maximum local 
Mach number, M|__^gj^, over the airfoil . (Saturation occurs when the 
pressure and temperature of the gas are on the vapor-pressure curve.) 
However, as reported in reference I, for the 0.3-m TCT, heterogeneous 
nucleation has only been detected below temperatures corresponding to 
freestream saturation. Therefore, if condensation effects occur over an 
airfoil in the 0.3-m TCT at temperatures above freestream saturati on, it 







complicated two-dimensional analysis. With regard to the applicability 
of the one-dimensional model, Wagner^^^ determined that the differences 
between one- and two-dimensional calculations for his condensation 
calculations are smaller than the differences between condensed and 
uncondensed flow. Wagner, therefore, concluded that the one-dimensional 
simulation should account for the main condensation effects in a 
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A 

and ^ Is the nondimenslonal area of the one-d1mens1onal nozzle with 
A* being the throat area at sonic speeds. By fol 1 owl ng an approach 
similar to earlier one-dimensional models of flow through a supersonic 
nozzle, as described In references 3, 4, and 5, the second step solves 
for condensatl on In the equivalent nozzle. In the present study, this 
step Involves Integrating the continuity, energy, and momentum equations 
along with a condensatl on equation based on the classical liquid droplet 
theory (CLDT) of homogeneous nucleatlon and Gyarmathy's droplet growth 
equation, using a fourth order Runge-Kutta 1 ntegratl on scheme. A1 so 
Incorporated Into this program Is the Beatti e-B r1 dgeman , ^ ® ^ equation of 
state which can be reduced to the Ideal gas equation of state. If 
desired, and a varl ety of suggested corrections to the CLDT - Tolman,^^^ 
nonisothermal Lothe and Pound, and that due to either Relss^^^^ 
or Kikuchl 

Data from a1 rfoll experiments In the Langl ey 0.3-m TCT and from 
earlier nozzle experiments are used to eval uate the various corrections 
to the CLDT, the real gas effects, and the sensitivity of the model to 
certain parameters. 
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EQUATIONS 
Flow Equations 

Since the flow over the supersonic section of an alrfolT can be 
approximated by the flow through an equivalent nozzle^ the condensing 
flow can be described following an approach similar to earlier models in 
references 3, 4, and 12. For this study, the flow is assumed to be one- 
dimensional, steady, frictionless at the boundary and adiabatlc—no heat 
transfer across the boundary. Instead of assuming an ideal gas as in 
earlier models, nitrogen is assumed to be a real gas which obeys the 
Beatti e-Bri dgeman ^ ® ^ equati on of state. The condensing droplets are 
assumed to occupy negligible volume with respect to the remaining vapor 
and have the same speed as surrounding gas. 

Continuity 

Since the flow is steady through the nozzle and the droplets are 
assumed to be at the same speed as the flow, the conti nuity equati on can 
be wri tten as 

Pj uA = m = constant (2) 

The total density, pj , is made up of two parts , such that 

Pf = Pq P (3) 

where p^ is the equivalent densi ty of the condensate dispersed 
throughout the same volume as the vapor density p . By solving for 
p 5 equation (3) can be expressed as 
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Py(l ” -^) - p (4) 

Since p' and pj are over the same volume, p' /pj can be replaced by the 
mass fraction of the condensate, g, which is def i ned as the ratio of 
the condensate mass to the total mass. Therefore , the total density 
pj can be expressed as 

By substituting equation (5) into equation (2), the continuity equation 
can be expressed as 


uA “ ft (6) 

The di f ferenti al form of the continuity equati on is then wri tten as 


p ix 



+ + 
u M 


1 dA _ 
A’Hx " 


0 


( 7 ) 


Momentum 

From the assumptions that the droplets move at the same speed as 
the gas, and that the flow is frictionless at the boundary of the 
nozzle, the momentum equation takes the form 


or 


u 






du 

W 


( 8 ) 
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Energy 

From the assumptions of no heat transfer across the nozzle boundary 
and steady one-dimensional flow, the energy equation can be obtai ned 
from the first law of thermodynamics as 

2 

f-+ h - g [Cp (T T^) + L(T^)] const (9) 

Because the first term in the brackets and the di f f erence between latent 
heat evaluated at T and are small compared to the rest of the 
equation, the energy equation can be simplified to the following form 
which is also used in references 2, 3, 4, and 12: 

u 2 

^ + h - gL * const (10) 

where L is calculated at the gas static temperature . The di ff erenti al 
form of the energy equation is then wri tten as 

,U, 

Equation of state 

To complete the solution of equations (7), (8), and (11), the 
equation of state is added to the flow equations. Because pressure is a 
function of density and temperature, the derivative can be expressed as 



( 12 ) 



The derivative of specific enthalpy of the gas, h, can be expressed as 


* . C 41 + (v - Kli) ) ^ 

dx p dx IT p OK 


(13) 


from reference 13. By substituting i for v and using the 
mathematical i denti ty 


(|E)^ (|p)p (|I)p = -1 


( 14 ) 


equation (14) can be wri tten as 


^ P P ' T ; 


(15) 


Solution of flow equations 

Because can be calculated from the nucleation and growth rate 
equations, to be di scussed later, and A is al ready determined, these 
values are treated as known quantities. Therefore, the 1 i near system of 
equations (7) , (8), (11), (12), and (15) remai n for the derivatives 

and . Because the 1 atent heat, L, is primarily 
a function of temperature, the derivative may be expressed as 


dL dL dT 

■ax " 


(16) 


By substituting equations (8), (15), and (16) into equation (11) and 
collecting terms, the energy equation can then be wri tten as 




i* "* y P \ Q I ^ p ^ P “I* 


/(H) ) +# (c 




dLv 


dq / • \ 
(17) 
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which can be expressed In the form 

By solving equation (12) for , substituting this result along with 
equation (8) Into (7), and collecting terms, the continuity equation can 
be written as 

.An r 1 ^ Iw dT r ^ (9p/3T)p__ _ 1 1 dft 


which can be expressed In the form 


■ 1. 

■1 “ ix 


+ E, 


dT 


= E, 


( 20 ) 


Using Kramer's Rule, equations (18) and (20) can be solved 
simultaneously for and In the following form 


du = i . ^3 h “ h h 
isix " u * 

dT ^ ^1 ^3 " ^3 ^1 
■3x Cj Eg - Cg' Ej 


(21) 

( 22 ) 


From the continuity equation (7), 


dp 

"Hx 


Is determined as 


dp 

•ax •" 



1 du . 1 dA% 

■— -y ' " » y -y*-} 

u ox A ^ 


( 23 ) 


Finally, Is calculated from the momentum equation as 



dg 

Hence, from equations (8), (21), (22), (23) , and calculated from 
the nucleation and growth rate equations to be discussed in the next 
section, the condensing flow over the airfoil can be solved. 


Homogeneous Nucleation 

The condensation over an ai rfoi 1 is assumed to be the result of 
droplets formed by the homogeneous nucleation of the nitrogen gas 
although the following equations are for any pure vapor. In the 
homogeneous nucleation process there is a size-dependent energy barrier 
which must be overcome such that all droplets with radius, r, less 
than some critical droplet radius, r*, will tend to evaporate while 
all droplets with radius larger than r* will be stable and grow. The 
critical radius is calculated from the Gibbs-Thomson^^^ or Kelvin^^^^ 
equati on 


£n ( 



2 0 


with the result that 


'■* " -pyTr-'a Wp--"!' 

where a is the surface tension of the droplet, is the droplet 
density, and Pgg^ is the saturated vapor pressure calculated at the 
temperature of the gas, T. At saturation p/Pgg^ “ 1 and r* » «; 
however, r* decreases as p/Pgat ^ increases and the probability 




increases that random collisions of gas molecules will result in a 
growing droplet with radius larger than r*. 

The formation rate of critical-sized droplets is represented by the 


classical liquid droplet theory, CLDT, which assumes that spherical , 
liquid droplets are formed which retain the bulk liquid properties and 



equation (26) can be rewritten as 


W = c/3 


(27) 


As menti oned in reference 12, this result is also true for the formati on 
of a solid cluster of molecules, or crystal of quite general shape 
where a is a "surface tension"-! ike term based on the free surface 


1 


since the Introduction of CLDT, many researchers have attempted to 
improve the theory by accounting for some of the physics that are 
disregarded by the basic CLDT. This secti on revi ews some of the best- 
known efforts to Improve the nucleatlon rate predicted by the classical 
theory. 
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To! man 

Because of possible error in assuming the planar value for surface 
tension of a liquid droplet, the effect of dropl et size on surface 
tension has been studied by several people.^^®^ One of the first to 
study the effect of droplet size on surface tension was Tolman^^^ who 
developed a correction to surface tension of the form 

0 " cr^/(l + ) (29) 

where a is the planar surface tension and 6 is a constant based on 
the intermolecular distances of the liquid which is on the order of 
several fluids. Tolman comments that less and less 
confidence can be placed in this expression as droplets get smaller and 
smaller. Typical , droplets for the present calculations contained 50 
molecules. Applying Tolman's correcti on to the surface tension, 
Sivier^®^ calculated the critical radius in equation (29) to be 

r* = r* (a ) - 26 (30) 

oo 

where is calculated from equation (24) using o^. Therefore, 

by using equation (29) and (30) and the appropriate Tolman constant, 6, 
the surface tension can be corrected and thus the CLOT can be modified. 

Lothe and Pound 

Because they felt that the free energy of formation in the CLOT, 
represented by the net work in equation (26), was not complete, Lothe 
and Pound^^^ proposed a correcti on to the CLOT which includes the free 




translational and rotational energies. Using the macroscopic or planar 

surface tension, a , the Lothe and Pound correction factor can be 
'00 

wri tten In the following form from reference 3: 
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“ ‘*’LP * ‘^CL 


(31) 


where 


2Trn* m k T 2/2 2 * 3/2 

♦lp«— ^ ] 


(32) 


T 


rep 


where 


n* = -j IT r*2 p^/m^ 


(33) 


and the moment of Inertia, I*, of critical spherical droplets Is given 
by 

I* = |- n* m^ r*2 (34) 


The repl acement factor, Q^gp, which compensates for the six degrees of 
freedom Introduced from the translational and rotational energies, 1s 
based on. the free energies of separation. For water, Qy.gp has been 
approximated between 10^ and 10® in reference 17. With Q^gp * 10^, 

(|>l_p has been estimated at 10^®. Because of the added 
translational and rotational energies, the Kelvin equation should be 
modified to include the correction term^^^^ 


iin (p/Pgg^) 


Z a 4 
" ■fiF’ 


(35) 



concluded that no rotational effects should be Included for a liquid 
droplet. In the appendum of Relss,^^^^ Reiss compares his and 
Kikuchl ' correction and concludes them to be essentially the same 
except that Kikuchl takes dropl et curvature dependence Into account. 
Disregarding curvature dependence, this correcti on will be labelled 
for el ther Reiss or Kikuchl and can be wri tten as 



so that the nucleatlon rate can be expressed as 

^ “ % '^CL 
Nonisothermal 

In contradiction to the assumption that the forming dropl ets are In 
thermal equilibrium with the condensing vapor which, as discussed by 






is valid only when a large quantity of carrier gas is present to 
di s si pate the heat of condensation, Wu, following Feder , et.al 
incorporated a correction factor to account for the temperature 
di f f erences between the droplets and the vapor that result from the 
necessary heat transfer between the f ormi ng dropl ets and the surrounding 
gas. When simplified for a pure vapor , this noni sothermal correction 
factor can be written as 

(38) 

T+1 + 2W - 1)[ - 7 - (1 - in 

SO that the nucleation rate can be expressed as 

J = *KI • JCL 

This noni sothermal multiplier reduces the nucleation rate by factors as 
large as 22 for the conditions analyzed herein. With regard to the 
growth equation proposed by Gyarmathy^^^^ discussed in the next secti on, 
a noni sothermal correction factor would not be proper because he assumes 
that the droplets and vapor are in thermal equilibrium at cri t1 cal 
conditions. 



that describes droplet growth in free-molecular, transition and 
continuum regimes is needed. As described in reference 22 , 


Gyarmathy developed a continuum equation which closely approximates 


droplet growth in the free-molecular regime as well . For a pure vapor, 
Gyarmathy' s equation takes the form of 



where is the Prandtl number of the gas, 5 is the thermal 
accommodation coefficient, which is assumed to be 1, and T^. is the 
temperature at the surface of the droplet expressed as ^ ^ 


T + (T 


sat 


T)(l - 


(42) 


where is the saturation temperature calculated at the pressure of 

the gas. The change in enthalpy, Ah^g , from the gas at T to the 
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liquid at is expressed as 



= Cp(T - T^) + 1. (y. 


( 43 ) 


By neglecting Cp(T - T^) and evaluating the 1 atent heat L at the 
temperature of the gas rather than in the energy equation (9) and 
in equation (43), the total temperature of condensation onset is 
typically reduced by less than 0.25K. Therefore, by approximating 
Ah. as L(T), equation (41) can be rewritten as 


dr _ X 
at pT 



,2/S^ 


1 (y7S~F P X Y+r^ 



(44) 


This dropl et growth equati on assumes no coagulati on of droplets. When 



Calculation of Condensate 

The condensate formed in the two-phase flow is determi ned from the 
integral of the condensate mass fraction g(x)^^^ , which is wri tten as 


g(x) = 1 p J(6)A(6) r^ (6.x) d6 (45) 

ifl -W 

where g, the nucleation rate J, and the effective nozzle area A are 
functions of position along the nozzle, x. The quantity r(6,x) is the 
radius at x of a droplet formed at 6, which can be written as 
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r(6, x) = r(6, S) + f - d e 

6 ^ ^ 


(46) 


where r(6 , 6) 1s the critical radius at 6 and where 


3r _ 1 dr 

W " uTT) lar 


(47) 


dr 


where Is determined from equation (44). Differentiation of equation 
(45) with respect to x gives 


® “ (~y pg J(x) A(x) r^ ( X , X ) + / 4irp^ J(6 ) A(6 ) r^(6 d6 ) 




(48) 


S f* 1 

where Is replaced by "if equation (47). Equation (48) 

shows that the Increase in mass fraction consists of the droplet mass 
fraction being created In dx and of the growth of particles created 
before the current x location. 

Using a fourth order Runge-Kutta Integration scheme , equations 
(48), (47), (21), (22), (23), and (8) are solved for the flow 
variables g, r, u, T, p, and p. The properties of gaseous and 
condensed nitrogen which are necessary In solving these flow equations 
are discussed In the next section. 
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PROPERTIES OF NITROGEN 

The properties of gaseous, liquid, and solid ni trogen are needed in 
predicting the formation of liquid or solid ni trogen droplets from 
nitrogen vapor. 

Gas Properties 

The properties of gaseous nitrogen are known relatively well down 
to the triple point (temperature, 63.148K, and pressure, .1237 atm, at 
which solid, liquid, and gaseous states of nitrogen coexist).^^^^ The 
extrapolated equations for the ni trogen gas properties can only be 
assumed to be valid for describing the nitrogen gas below the triple 
point and the metastable gas which exi sts in the liquid regime before 
the gas condenses. 

Equation of state 

The Beattie-Bridgeman (B-B) equation of state for ni trogen is used 
herein to i ntroduce real gas effects. This equation of state is 
surprisingly accurate and is represented by^®^ 

p = pRTd - ^ p)[l + Bq p(1 - bp)] - AqP^(1 - ap) (49) 

where the values of the constants for nitrogen are listed in table 1. 

By setting these constants to zero, the ideal gas equati on of state is 
obtained 

p » pRT (50) 

Using the B-B or ideal gas equation of state, the following 
expressions are used for calculating the specific enthalpy, h, 
specific entropy , s, heat capacities, Cy and Cp, and sound speed. 
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a, for nitrogen gas 


.( 6 ) 


h-h*T ll c“ dT (51) 


0 T 0 p p 
0 


P T 


^ dT - R in (^) +jj^ 

p 


S " S° "T^o h 


(52) 


(53) 


c = c, . T (|^)^ /Cp2(M),] 


(54) 


a = [J(|£) 

Cy 3P I 


(55) 


where h*j - 309500 J/kg and s*t = 6836.09 J/kg-k are reference 
'o 0 

values of specific enthal py and entropy, the reference pressure, pQ, 

is 101325 N/m^, and the reference temperature, Tq, is 298.15K.^^^^ The 
zero pressure values of specific heats, Cp® and Cy® are represented by 
the following ideal expressions 


®n° = R 
P r-^ 


(56) 


where the ratio of specific heats , y = -r^, for "ideal “ ni trogen or any 

^v 

ideal diatamic gas, is 1.4. Therefore 


0 = 


3.5R 


(56a) 


and 


r 0 = r 0 _ R 
Cy Cp K 


( 57 ) 
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In figure 3, the results of calculating Cp using the B-B equatl on 
of state (eq. (50)) are compared with the constant Ideal gas value. 

A1 though not shown, similar resul ts would be calculated for Cy and 
y. As in the fol 1 owl ng f1 gure the gas properties are calculated along 
the vapor pressure, V-P, curve (saturation boundary) and constant 
pressure lines of 0.1 and 2 atms. The constant pressure lines, which 
are the range of local static pressures typically encountered during 
airfoil tests, extend from the gas regime across the V-P curve and into 
the metastable gas regime by about 15K, which would be below the 
temperatures at which the onset of condensation was observed during the 
airfoil experiments to be discussed In a later section. 

Mean free path 

The mean free path, I , Is the average di stance travel ed by a 
molecule between collisions. The ideal expression assuming a rigid 
spherical molelcule^^^^ Is used for I as shown 


vT TT 



(58) 


where ~ Is the molecular volume and d is the rigid sphere d1 ameter . 

By substituting d = 3.75 x 10“^®m for nitrogen,^^®^ the mean free path 
for gaseous nitrogen can be wri tten as 


- . 7.44 X 10”^ 


P 


(59) 
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Thermal conductivity 

The thermal conductivity of nitrogen, X, can be calculated using 
equation from Jacobsen' 'which is made up of two terms 

X = X^(T) + A X^(p) (60) 

where the dilute gas contribution, x^, and the excess, or dense, fluid 
contribution, ax^, which is basically a function of density only, are 
1 i sted in table 2. 

In Jacobsen,^^^^ there is an additional term which is the 
enhancement due to the influence of the critical point. For the present 
study, this term is del eted because it is generally small and because 
the high pressures (~ 30 atms) in the neighborhood of the cri ti cal 
poi nt are not typical of condensation onset in transonic cryogenic wind 
tunnels. Equation (60) is compared over the temperature range of 
interest with a Sutherland-type expression used by Wagner^ which is 
given as 

X = 2.067 X 10"^ j (61) 

As shown in figure 4, the results of both equations are similar but the 
X calculated with equation (60) is slightly pressure dependent. For 
the present studies, the Sutherland-type expression is used. 

Viscosity 

The viscosity, n , of nitrogen is calculated from an equation in 
Jacobsen, which contains two terms and is wri tten as 
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n = n^d) + A Ti^(p) (62) 

v/here the dilute gas contribution, n., and the excess, or dense, fluid 
contribution, Hg, which is primarily dependent on density, are 1 i sted 
in table 3. In Jacobsen, there is another expression for An^ for 
when the gas density p is greater than 800 kg/m^, but this term is not 
included because gas densities of this magnitude are not typical of 
condensation onset in transonic cryogenic wind tunnels. Equation (62) 
is compared over the temperature range of interest with a Sutherland- 
type expression used by Wagner, which is given as 

n * 1.378 X 10“® (63) 

For the range of local pressures during the airfoil tests, the results 
of both equations are similar for the temperatures above 70 K as shown 
in figure 5, but the n calculated from equati on (62) is slightly 
pressure dependent. At the lower temperatures, the resul ts of equation 
(62) are questionable and, in fact, the Sutherland-type expression for 
viscosity is recommended in place of the Jacobsen expression. 

Prandtl number 

The Prandtl number, Pr, is defined as 
Cj, n 

Pr = (64) 

For constant pressures and the 'Jacobsen expression for thermal 
conductivity, X (eq. (60)), and viscosity, n (eq. ( 62 ) ) , and the 
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Ideal value of Cp = 3.5R, the Prandtl number approaches 0.73 for 
temperatures above 70K. This value agrees with the constant value for 
Pr of 0.72 (see fig. 6a) calculated using the Sutherland- type 
expressions for x (eq. (61)) and ri (eq. (63) ) from Wagner^ and 
assuming the Ideal value for Cp of 3.5R. However, at the lower 
temperatures, Pr 1 ncreases significantly from the Sutherland value 
of 0.72, because of the increase In Jacobsen's expression for viscosity 
at the lower temperatures (see fig. 5). 

Using the B-B equation of state and thus the real -gas value of 
c„, the Prandtl numbers calculated using the Jacobsen and Sutherland 

H 

expressions for x and n are in good agreement with each other at 
temperatures above 70K, where both calculations are greater than the 
1 deal gas value of Pr (see figs. 6b and 6c) because of the real gas 
value of Cp which Is shown In figure 3. Even though, at the lower 
temperatures, the real gas values of Cp approach the 1 deal gas Cp of 
3.5R, the Prandtl numbers calculated from the Jacobsen and Sutherland- 
type expressions for X and n again do not agree because of the low- 
temperature differences In r\ (see fig. 6b and 6c). The Prandtl number 
calculated from the Sutherland-type equations for X and n might be 
appropriate for both above and below the triple point because of (1) the 
good agreement between the Jacobsen and Sutherland- type expressions for 

n above 7 OK and X over the full temperature range for constant 
pressures and (2) the more reasonable n values calculated from 
Sutherland-type expression at temperatures below 7 OK although there Is 
no d1 rect experimental evidence at temperatures below the triple pol nt . 
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By equating equations (70) and (71), solving for L, and replacing v 

1 

with — , the latent heat of sublimation can be expressed as 

P 


b p 


sat 


-) tn 10 


(72) 


where Pgg^ and b are from reference 29. As can be seen in figure 9, 
a relatively constant value of 2.43x10'' J/kg for the nitrogen 1 atent 
heat of subl imation is calculated from equation (72) for both ideal and 
real gas equation for p . This value agrees well with Wu's^^^ value of 
2.42x10^ J/kg for the ni trogen latent heat of sublimation. Griffin^^^ 
and Sivier extended the latent heat of vaporization below the triple 
point. Extension of equation (69) , the curve fit to Jacobsen's values 
of latent heat for the liquid, below the triple point in figure 9 shows 
that a markedly different value of nitrogen latent heat of sublimation 
could result. Therefore, the current program uses a constant value of 
2.43x10® J/kg for the nitrogen latent heat of sublimation. 

Surface tension 

Surface tension is a famil iar concept for liquids and can be 
readily measured. For 1 iquid nitrogen, the surface tension can be 
expressed as^^^^ 


a = 0.0297074 (1 - T^ /126.2)^‘^^^^^ 
while earlier expressions such as the one by Sivier^®^ 


(73) 


0 = .02394 - .0001933 T„ 


(74) 















SI 


3 2 

of a /p^ is an important quantity to predicting the nucleation rate. 
Consequently, knowing the values of the surface tension and density used 
in a particular paper is necessary to properly compare the di f f erent 
theoretical expressions different researchers have used for predicting 
the onset of homogeneous nucleation. In other words, many differences 
in theoretical expressions may be overshadowed by variations in 








local Mach number, I^L-max* over the airfoil was approximately 1.4 
for the total pressure range tested. 


Compari son of the CAST- 10 airfoil data with the one-dimensional 


computer model discussed in the previous sections is restricted to the 


supersonic region over the airfoil , because this is where condensation 


due to homogeneous nucleation will occur first (see fig. 12) and because 
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Dankert 

The most useful data of Dankert^ for the present comparl sons are 
his free- jet nozzle data from his figure 34 in which p^ = 3.0 atm and 
the throat di ameter is 0.005m. For = 125K, 160K, and 17 5K, the 
computer model dramatically underpredicts condensation effects using 
CLOT, CLDT-T, or CLDT-DK which are not shown in figure 18. In figure 
19, the computer model using CLDT-LP with Q^.gp «= 10^ agrees well with 
data for = 125K and reasonably well with data for » 160 K and 
175K before a significant amount of condensati on occurs, which causes 
the assumptions used in the model to become invalid. In f i gure 18b, 
onset of condensation effects are overpredicted for all three cases 
using CLDT-WU except when large amounts of condensate are present, where 
the assumptions used in the model are inaccurate. In the next section, 
further compari sons are made with wedge nozzle data. 

Nagamatsu and Willmarth 

Because Sivier's computer model in reference 5 agrees with the 
wedge nozzle data from figure 14 of Nagamatsu and Willmarth^^^^ for run 
9-5 with p^ » 8.21 atm and = 295K, the present computer model Is 
compared with these data. No condensation effects are predicted using 
CLOT, CLDT-T, CLDT-LP , CLDT-DK, or CLDT-SIV. Condensation effects are 
also underpredicted using CLDT-WU as shown in figure 19. This disagree- 
ment between the computer model and the Nagamatsu and Willmarth data 
could he due to either impurities in the nitrogen used by Nagamatsu and 
Wilmarth or in the assumptions in the computer model . In particular, 
since Sivier's properties for the solid nitrogen density and surface 
tension are used in CLDT-SIV, the primary di f ference between Sivier's 
model and CLDT-SIV developed herein is that Sivier assumes the droplets 





model . Homogeneous nucleatlon is represented by the CLDT-RK for liquid 
nitrogen and the CLDT-LP with Q^gp = 10^ for solid nitrogen, although 
more experimental comparison data are needed to verify these choices. 
Growth of the droplets is represented by Gyarmathy radius- 
dependent droplet growth equation. In fi gure 20, the predicted onset of 
condensation effects is said to occur when the pressure deviates from 
the no-condensation pressure by 0.25 percent at x/c = 0.25. This is an 
arbitrary location ahead of the recompression shock for the CAST-10 
experiment, which is the pressure distribution used in predi cti ng this 
onset curve. For other pressure distributions, a d1 f f erent 1 ocati on may 


P 
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'liquid nitrogen condensation agrees well with the experimental onset 
which was measured from the CAST- 10 data. Bel ow = 3 atm, exact 
prediction of homogeneous nucleation may be unimportant for the 0.152m 
CAST-10 airfoil because the predicted onset of condensation falls below 
freestream saturation, where condensation on pre-existing seed particles 
(heterogeneous nucleation) can occur upstream of the airfoil and 
influence the aerodynamic data at temperatures above which homogeneous 
nucleation could occur. The inflection of the predicted onset curve 
at p^ = 2.5 atm is due to part of the flow being above and part of the 
fl ow being bel ow the triple point during the calculation. 



underpredicts condensation effects. 

5. With the Beattle-Bridgeman equation of state used In calculating 
the expansion Isentrope, the change In condensation effects Is 
relatively small. 

6. The calculations are very sensitive to the value of surface 
tension. 
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7. With the Reiss or Kikuchl correction to the classical theory, the 
computer model Indicates that below = 3 atm for the 0.152m 
CAST-10 airfoil with = 0.65 and a = 6°, exact prediction of 
homogeneous nucleatlon of nitrogen may be unimportant because the 
predicted onset of condensation falls below freestream saturation, 
where heterogeneous nucleatl on can occur upstream of the airfoil 
and Influence the aerodynamic data at temperatures above which 
homogeneous nucleatlon could occur. 

8. To further val 1 date the onset curve, more data need to be taken. 
These data must be taken very acurately because of the high 
sensitivity of the onset of condensation to total temperature. 

In predicting the conditions which would give rise to effects due 
to the condensation of solid nitrogen, the results of comparl ng the 
computer model with experimental data are summarized In the following: 

1. The surface tension and density are both Important when comparing 

the nucleatlon rates of solid nitrogen because of the dependence of 

3 2 

the nucleatlon equation on o /p . 

2. Using Wu's and Sivler's expressions for solid nitrogen surface 
tension and density, the classical theory dramatically overpredicts 
condensation effects for the CAST- 10 airfoil data. 

3. A theoretical expression for the surface tension of solid n1 trogen 
Is calculated from a modified version of Tabor's expression for 
surface tension. 
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4, Use of the classical theory with the Lothe and Pound correction 

Q^gp = 10^, density from the Scott data and the expression for 
surface tension mentioned above, agrees with the Dankert free- jet 
data while use of the classical theory with Wu's surface tension 
and density overpredicts the onset of condensation effects. 

5. Comparison of the present computer model with the Nagamatsu and 
Wlllmarth nozzle data and Sivler's computer model Indicates that 
assuming droplet temperatures equal to satuatlon temperature rather 
than Initially at gas temperature could greatly Increase nucleatlon 
rate, although this assumption would be controversial . 

6 Before any definitive conclusions can be made concerning onset of 
solid nitrogen formation, further studies on the properties of 
solid nitrogen and more experimental data are needed. 
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APPENDIX 


MODIFIED TABOR THEORY FOR SURFACE ENERGY OF SOLID NITROGEN 
Tabor^^^^ outlined an approach which related the 1 atent heat Of 
sublimation to the free surface energy. Basing his expression on the 
[100] face of a face-centered cubic (fee) solid, he deri ved the 
following equation for surface energy 


.267 (— |-) ( 


'■a ^s, 


2/3 


A-1 


where the constant .267 changes by small amounts for other faces and 
structures.^^®^ Therefore, because 3 - phase solid nitrogen is a 
hexagonal closed-packed (hep) solid, equati on A-1 must be modified. 
Assuming that the dropl ets form along the lowest energy surface which is 
the [0001] face for a hep solid, the surface energy is recalculated 
following the basic approach outlined in Tabor 

On the [0001] face of a hep solid, each surface molecule has 9 
"nearest neighbors", whereas inside each crystal structure, each 
molecule has 12 nearest neighbors. If each crystal contains a total 
of N atoms with Ng surface atoms , the total bond energy of the 
crystal , counting only nearest-neighbor interaction is 


j Ac [12(N-Ng) + 9 Ng] = I Ac [12N - SN^] 


A-2 


where Ac is the potential energy between each molecule and its 
neighbor, and where the factor of ^ is introduced to avoid counting 
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1 

If there were no surface atoms, the energy would be ^ Ae 12N. 

3 

Therefore, the energy is reduced by j Ae Because Ae is 

negative, this is a positive increase in energy which can be associated 
with the surface energy. If the surface area of the crystal is A, the 
surface energy can be identified as 

oh “ ^ A e Ng A“3 


Therefore, equation A- 3 can be solved for o 



s _ 3 

~X " 7 


Ae z 


^-4 


where z is the surface molecules per m^. 

The latent heat of sublimation can be expressed in terms of Ae • 
for a closed pack structure, such as a hep structure, by the expression 




(Note: In Tabor's method L is in J/mole, but for equation A-5, 

is di vi ded by the molecular weight, p , because L has units of J/kg.) 
Equation A-5 can be solved for Ae as 

L p 

Ae = A-6 
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and then Ae can be replaced in equation A-4 so that 


Lg P 


'■s “ , 



center of its nearest neighbor. If y is the molecular weight and pg 
the density, the volume of each molecule is also represented by 



Compari son of equation A~10 and A-1 shows that the constant .229 in 
A-10 is approximately 15-percent lower than the constant in A-1. 
Therefore, equation A-10 would predict a lower surface energy than 
Tabor's expression. 
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TABLE 1.- BEATTTE-BRIDGEMAN EQUATION OF STATE CONSTANTS FOR NITROGEN, 

REF. 6 

Aq = 173.60 N-m^/kg2 
a = .0009342 m^/kg 
Bq = .001801 m^/kg 
b = - .000247 m^/kg 
C = 1499. m^k^/kg 



J 

TABLE 2.- THERMAL CONDUCTIVITY FOR NITROGEN, "i-K-sic , 

JACOBSEN REF. 24 

X = Xq(T) + A Xpjp) 

X = I t(^"3) 

° 1=1 

AX^ = .00002195902219 p + .006387370699 ^ - 1.0) 

Ai = -6.8939127475 
A 2 = 3.5226118983x10“^ 

A 3 = -6.8357539823 x lO"^ 

A 4 = 1.5832717315 x 10“^ 

A 5 = -2.6418423047 x 10"^ 

Ag = 3.6093309138 x 10“1° 

A 7 = -2.5555598476 x lO’^^ 

Aq = 8.5635041641 x 10'^^ 

Ag = -1.0717599406 x 10“^0 
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TABLE 3.- VISCOSITY FOR NITROGEN, 
JACOBSEN REF. 24 


N-sec 



» 


n = (p) 

„ C, T 

0 i=l ^ 

Cl = 7.4165322904 x 10'^ 

C 2 = -1.5834400475 x 10"^ 

C 3 = 3.8530771011 x lO''^ 

C 4 = 8.0133713668 x 10"® 

C 5 = -8.9203123846 x 10"^^ 

Cg = 8.9059711315 x 10"^^ 

C 7 = -5.3779372664 x lO’^^ 

Cg = 1.7398277309 x lO'^O 

Cg = -2.3084044942 x lO"^^ 

7 .. 

Ann = I D. (O.OOlp)^ 
i=l ^ 

Di = 2.3083514362 x 10"® 
D 2 = -9.3636207171 x 10”® 
D 3 = 9.0339186452 x 10”^ 
D 4 = -4.1832067163 x 10”® 
D 5 = 1.0897627893 x 10’® 
Dg = -1.2913856376 x 10’® 
D 7 = 5.9782049913 x 10"® 



Reference 

Range 

a 

b 




TABLE 4.» VAPOR PRESSURE CURVE 


^at = 


(a-b/T) 


Dodge and Davis Ref. 29 
63.148<T<126.20 K 
3.93352 
304.494 sec 
101325 N/m^ 


Frels, et. al . Ref. 30 
35.5<T<63.148 K 
7.614676 
356.281 sec 
133.32 N/m^ 


c 
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TABLE 5.- COMPARISON OF PREDICTED CONDENSATION ONSET TOTAL TEMPERATURE 
CAST-10 AIRFOIL. M = 0.65, p. = 5.0 atm 

CO 

BASELINE-CLDT, T^ = 99K. * 0.0266, x/c - 0.25 


Modification 

Onset T^, K 

T^ " 99K 

'^/'^CL 

CLDT 

99 

0 

1 

CLDT-T, 6 = 0.25xl0"^°m 

100.7 

1.7 

500 

CLDT-RK 

100.3 

1.3 

no 

CLDT-LP, Q^ep = 10^ 

103.8 

4 «8 

loii 

Nonisothermal 

98 

-1.0 

0 

CLDT, a reduced 10% 

102 

3 

10® 

CLDT, B-B equation of state 

, 99.8 

.8 

18 




AIRFOIL 



Figure 2.- Approximating airfoil pressure distribution by 
an equivalent one- dimensional nozzle. 
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Figure 3.- Comparison of specific heat at constant pressure as a function of temperature 
plotted for constant pressures of 0.1 and 2 atm and for the pressures associated with the 
vapor- pressure curve using the Beattie- Bridgeman equation of state with the ideal gas value. 
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Figure 4 .- Thermal conductivity of nitrogen gas. 
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Figure 5.- Viscosity of nitrogen gas. 



a) 


Both Jaccfcsen and Sutherland- type \ and n and ideal gas c 
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Figure 6.- Prandtl 
confined with c 

P 


nuntier calculated using either Jacd>sen or Sutherland- type A and n 
from either Beattie- Bridgeman or ideal gas equation of state. 
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ro 
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Figure 7.- Vapor- pressure curve for nitrogen. 
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Figure 8.- Condensate density for nitrogen. 
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VAPORIZATION 



Figure 9.- Latent heat values for nitrogen. Line labelled IDEAL CALC, 
uses ideal gas equation of state values for gas density while line 
labelled REAL CALC, uses Beattie- Bridgeman equation of state values of 
gas density. 
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Figure 10.- Surface tension for nitrogen. 
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Figure 11.- Comparison of different expressions for a /p^. 




Figure 12.- Pressure coefficients in region of interest for 

CAST- 10 airfoil. M = 0.65, a= 6°, and c= 0.152m. 

00 ’ 
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T^. K 

Figure 13.- Saturation lines as well as lines showing total conditions 
for CAST- 10 airfoil at a=60. O -Experimental onset of condensation 
effects. 
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Figure 15.- Sensitivity of calculations to small 
changes using CLDT-RK, : 
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Figure 17.- Comparison of CAST- 10 datasO, and 
theory for R^= 15x10°, M^= 0.65, and a- 6°. 




a) Comparison with CLDT-LP, Q 
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pure 18.- Comparison of Dankert data and theory. Pressure ratio 
as a function of orifice diameters downstream of orifice. 


p.= 3.0 atm, D= 0.005m. 



b) Comparison with CLDT-WU. 


Figure 18.- Continued. 



x/l 


Figure 19.- Compari son of Nagamatsu and Willmarth data and 
CLDT-WU. Pressure ratio as a function of position. 
p^= 8.21 atm, 295 K, and 1= 0.0254m. 
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Figure 20.- Predicted onset curve for 0.152m CAST- 10 airfoil. 0.65, 
and a= 6°. O -Experimental onset of condensation effects. 
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